arXiv:1501.01952vl [cond-mat.mes-hall] 8 Jan 2015 


Complete All-Optical Quantum Control of Electron Spins in InAs/GaAs Quantum 

Dot Molecule 

Guy Z. Cohen 

Department of Physics, University of California, San Diego, La Jolla, California 92093-0319 

(Dated: January 9, 2015) 

The spin states of electrons and holes confined in InAs quantum dot molecules have recently come 
to fore as a promising system for the storage or manipulation of quantum information. We describe 
here a feasible scheme for complete quantum optical control of two electron spin qubits in two 
vertically-stacked singly-charged InAs quantum dots coupled by coherent electron tunneling. With 
an applied magnetic field transverse to the growth direction, we construct a universal set of gates 
that corresponds to the possible Raman transitions between the spin states. We detail the procedure 
to decompose a given two-qubit unitary operation, so as to realize it with a successive application 
of up to 8 of these gates. We give the pulse shapes for the laser pulses used to implement this 
universal set of gates and demonstrate the realization of the two-qubit quantum Fourier transform 
with fidelity of 0.881 and duration of 414 ps. Our proposal therefore offers an accessible path 
to universal computation in quantum dot molecules and points to the advantages of using pulse 
shaping in coherent manipulation of optically active quantum dots to mitigate the negative effects 
of unintended dynamics and spontaneous emission. 


I. INTRODUCTION 


Spins in semiconductor quantum dotsEla (QD) are 
promising candidates for satisfying- the DiVincenzo 
criteria for i-quantum computations such as- state 
initialization, m3 coherent spin manipulation,and 
qubit-specific measurements.EJJO The advantages of QD 
spin qubits over other qubit realizations are numerous. 
Quantum dots are compatible with existing semiconduc¬ 
tor technology and can be grown in millions in ..a reg¬ 
ular 2D array on a single millimeter-sized chip£3 QD 
spin qubits can_bc_integrated on a chip with photonic 
crystal cavities, E jlj have shpH optical recombination 
and photon emissioH-,times,E30 can—be-, manipulated by 
fast single-qubi1nli3’t®l and two-qubitc3’E3 quantum gates, 
and can be entangled with adjacent qubits by tunnel¬ 
ing interactioncd or with, remote ones via entanglement 
swapping with photons.t5c3 Recently, great progress was 
made in developing methods to increase-the QD spin de¬ 
coherence time either through spin ecliocJ or suppression 
of nuclear-spin fluctuations.EH 

The study of the interaction mechanisms between QDs 
can lead to finding new ways to successfully manipu¬ 
late their quantum state, an important goal of quan¬ 
tum information technologies, like quantum computing 
and quantum cryptography. Such coupling can,-be ob¬ 
tained in optically active self-assembled QDs, I 
techniques epating for ultrafast laser i 
measuremeniGrEJ and coherent manipulationEKEddiU of 
the spin qubits in these QDs. Many coupling mecha¬ 
nisms of optically-active QDs were studied including elec¬ 
tron tunneliagjH3’E3 hole tunneling, ErEil exciton-mediated. 
interactional and electron-hole exchange interaction.!^ 
In particular, it was found that embedding two coupled 
QDs, which form a QD molecule (QDM), in a Schottky 
diode structure enablesrthe tuning of the relative energy 
levels of the two dots.E3E3 Like real molecules, QDMs can 



display bonding and anti-bonding states, that are sym- 
metriC|_and anti-symmetric superpositions of localized 
states.Ea However, unlike natural diatomic molecules, 
which always have a bonding ground state, QDMs can 
be tailored to have-molecular ground states with anti¬ 
bonding character .til 


With the QDM embedded in a Schottky diode, the 
tuning of the voltage bias .determines the stable ground 


portant is the 
where each of 


state charge configuration.GI Es 
case of a doubly charged QDMjI 
the QDs is charged with a single electron or hole, since 
then the QDM can be used as a double spin qubit system. 
An ambitious goal in quantum computation is complete 
quantum control of such a multiple-qubit system. This 
control can be achieved by a universal set of quantum 
gates, i.e. a set of gates with which an arbitrary uni¬ 
tary operation can be constructed. Many universal sets 
of quantum gates were proposed. Ref. ||| showed the set 
of all two-level unitary gates is universal. Then, Ref. |39| 
showed this set of gates can be implemented by combina¬ 
tions of single-qubit gates and the two-qubit controlled- 
NOT (CNOT) gates. Other universaLsets consisting of 
only two-qubit gates were also found .Grcll 


In this work, we present a scheme for complete quan¬ 
tum control of a two-qubit system, a QDM composed 
of two vertically-stacked singly-charged InAs QDs sepa¬ 
rated by a GaAs/AlGaAs tunnel barrier and embedded 
in a Schottky diode. The two electrons interact by ki¬ 
netic exchange, which is based on coherent electron tun¬ 
neling and evidenced by a singlet-triplet splitting of the 
QDM energy levels. Quantum optical control is realized 
through real or virtual excitation of the bottom QD elec¬ 
tron by laser pulses. The universal set of gates consists 
of the 5 two-level unitary operations that correspond to 
the possible Raman transitions in the Voigt geometry. 
We find that every two-qubit unitary operation can be 
realized through at most 8 such Raman transitions, and 
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that pulse shaping the laser pulses can substantially in¬ 
crease the fidelity of the operation. As a case in point, 
we show a reasonably high fidelity realization of the two- 
qubit quantum Fourier transform (QFT) in this system 
taking light hole mixing, decay and decoherence into ac¬ 
count. The results indicate the promise of QDMs as a 
platform for quantum computation and the potential of 
pulse shaping for increasing the fidelity of quantum state 
manipulation in optically active QDs. 

The paper is organized as follows. In Sec. || we present 
the system model and Hamiltonian. We write the Hamil¬ 
tonian relative to an orthogonal basis, find the eigen¬ 
states and eigenenergies, and plot the allowed transitions 
in the Faraday and Voigt geometries. In Sec. 0 we de¬ 
scribe a universal set of gates for the QDM. We give an 
algorithm for decomposing a given arbitrary two-qubit 
unitary operation to a product of operations realized with 
these gates and give criteria for a subset of these gates 
to still be universal. In Sec. 0 we present the Hamilto¬ 
nian of the interaction with the electromagnetic field and 
discuss the pulse shapes used to implement Raman tran¬ 
sitions and the effects of light hole mixing. In Sec. [v| we 
consider the effect of decay and decoherence. We write 
the Lindblad equation for the system and derive the ex¬ 
pression for the fidelity of an operation relative to an ideal 
one. In Sec. VI we show, for realistic parameters, how to 
apply the scheme to realize the two-q ubit QFT with rea¬ 
sonably high fidelity. Finally, in Sec. VII we discuss the 
key results and consider directions for future research. 




FIG. 1. Schematic diagram of the band structure of two ver¬ 
tically stacked self-assembled quantum dots embedded in a 
Schottky diode. The growth direction is the z direction. The 
height of the bottom/top dot is /is// it, while the height of 
the interdot tunnel barrier is do- Optical control of the two- 
electron spin state is achieved through real or virtual transi¬ 
tions to the optically excited states using laser fields tuned to 
create an exciton at the bottom dot only. The electric field 
in the a direction, F, is exaggerated for clarity. 


II. SYSTEM MODEL AND HAMILTONIAN 


The system under consideration consists of two InAs 
vertically stacked self-assembled quantum dots separated 
by a GaAs/AlGaAs tunnel barrier and embedded in a 
Schottky diode. The voltage of the diode is adjusted to 
the charge stability region where each dot occupied with 
a single electron at the ground state. The coherent tun¬ 
nel coupling of the electrons in the two dots is manifested 
in inter-electron kinetic exchange interaction, which gives 
rise to singlet and triplet electron states delocalized over 
both dots. The quantum states of the electrons are op¬ 
tically manipulated via laser fields that generate real or 
virtual electron-hole pairs (excitons) at the bottom dot. 
The magnetic field is initially taken as zero. 

With the growth direction taken to be in the z direc¬ 
tion, the potential experienced by a single electron in the 
conduction band is given by 

U(r) = ^mj_u 2 ± (x 2 + y 2 ) + U(z), (1) 

where m± is the effective electron mass in the x and y 
directions, approximated as independent of z, to± is the 
effective frequency of the parabolic confinement in the 
x and y directions, and U(z) is plotted in Fig. [j]. The 


Hamiltonian for the ith electron, in turn, reads 


hi — 


p? 


- p 2 


pL 


2m± 


2 m(Zi) 


U(n), 


( 2 ) 


where and p 4 (i = 1,2) are the position and momen¬ 
tum of the ith electron, and where the effective mass in 
the 2 direction is m(z). The two-electron Hamiltonian 
is constructed by writing Eq. (|j) for each electron and 
adding the Coulomb interaction term. It is 

'LL = hi + /i 2 + -j———r, (3) 

|ri - r 2 | 

where e$ = e 2 /(47re 0 )/{[e(2:i) + e(z 2 )]/2}, e being the 
electron charge, eo the vacuum permittivity and e(z) the 
relative dielectric constant. 

We consider the two dots to be separated such as to jus¬ 
tify a variational treatment in terms of atomic-like single¬ 
particle states localized at the individual dots. These 
single particle states are |H)| t)> 1-6)1 I), |T)| t) and 
|T)| j.), where \B)/\T) is a state corresponding to the so¬ 
lution of the Schrodinger equation for the single-particle 
Hamiltonian in Eq. (|2j) with U(z) containing only the 
bottom/top dot potential well. Listing the two-electron 
combinations of these states gives the product state ba¬ 
sis: | t,t), I t4), I 4-,t), I 14), I 4t,0) and |0,|t), where 
IcsjCTt) is the state with the bottom/top dot occupied 






















3 


by an electron with spin projection ctb/ot in the z di¬ 
rection, 0 denotes an unoccupied dot, and |t denotes a 
doubly-occupied dot. 


new basis reads 


+} = 2 -1 / 2 (| 0) + 0,tt))> 

(4) 

\S) = 2 -1 / 2 (| til) — l)t»> 

(5) 

|T 0 ) = 2- 1 /2(| t,|) + U,t)), 

(6) 

|T+) = lt,t), 

(7) 

|r-) = ll,l), 

(8) 

|-) = 2" 1 / 2 (||t,0)-|0,|t)), 

(9) 


With foresight, we choose a new basis composed of 
linear combinations of the product state basis states. The 

I 


where |±) are the symmetric/antisymmetric doubly oc¬ 
cupied states, IS 1 ) and \T 0 ) are the singlet and triplet 
Heitler-London states, and |T±) are the two other triplet 
states. The spatial parts of the triplet states are anti¬ 
symmetric, while those parts in the rest of the states are 
symmetric. This corresponds to the former/latter being 
anti-bonding/bonding states. 

The Hamiltonian in Eq. (||) in the basis of Eqs. (f|j9j) 
is 
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( 10 ) 


where I is the identity matrix, and the following defini¬ 


tions are employed: 

€B = (B\h\B), (11) 

e T = {T\h\T), (12) 

t = -(B\h\T), (13) 

tB = t — ubbbt — {B\T)eB, (14) 

tr = t — utttb — {B\T)ct, (15) 

J = ubtbt - 2(B\T)t. (16) 


The single-particle Hamiltonian h in Eqs. 0® is given 
in Eq . (p|). The matrix elements of this Hamiltonian in 
Eqs. ( |11| - |12| ) are the single-particle energies for the states 
| B) and |T), apart from small contributions from the 
other quantum well in U(z). The matrix element t in 
Eq. (Il3j) is termed the tunneling matrix element. The 
Coulomb integral Uijki ( i,j , k,l = B or T) is given by 

Uijki = f dndr 2 -, ———r^*( r i)<( , i( r 2)0/c(r2)^( r i)) 

J | ri — r 2 | J 

(17) 

with cj>i (r) being the atomic-like single-particle wave func¬ 
tion of the state |i). For the direct Coulomb integrals we 
use the abbreviations Uij = Uijji and Ui = uuu. 

We make the basis of states in Eqs. (Q-^[) orthogonal 
by replacing |+) and \S) in Eqs. (f^|]) by (| S) ± |+)). 
We then normalize the states of the orthogonal basis and 


write the Hamiltonian in Eq. © in the orthonormal 
basis. The diagonalization of this Hamiltonian shows 
the three triplet states \Tq), \T±) are degenerate eigen¬ 
states with an energy independent of the electric field 
F. The singlet state | S) is shifted in energy by the ki¬ 
netic exchange interaction below the triplet states, and 
this interaction also results in an admixture of the dou¬ 
bly occupied states |±) in the eigenstate dominated by 
| S). For the specific experimental values of hs =2.6 nm, 
hr = 3.2 nm, do = 9 nm (see Fig. [lj) and As_t = 125 
/xeV, taken from Ref. |2^, we plot the eigenenergies of Ti. 
as a function of the external electric field F in Fig. [|. 
The working point is taken as the one for which the 
singlet-triplet splitting is As-t = 125 /xeV, and at this 
point the singlet eigenstate is numerically obtained as 
\S) = 0.973|S)-0.173|+)+ 0.149|-). We term this state 
together with the triplet states, the spin states. 

We now consider the optically excited states in our sys¬ 
tem, namely the X 2 ~ states which consist of a negative 
trion at the bottom dot and an unpaired electron at the 
top dot. These states are | 4.t"fht)> I Itthl}) I lt-D->t}> 
and | lt-0-,1}, with | 4_tlt) = 2" 1 / 2 (| fi) - | It))I l!) 
and | |tl) = 2-V 2 (| ft) - I It)) I I) representing tri- 
ons, and with | j)) = ||, |) and | l)) = ||,—|) being the 
heavy hole states with 3/2 and —3/2 spin projections 
along z. The spin states can be excited to the optically 
excited states, which we term the trion states, through 
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FIG. 2. Calculated spin state energy levels for the quantum 
dot molecule system vs. F, the applied electric field in 2 di¬ 
rection, with the experimental parameters of of hs = 2.6 nm, 
hr = 3.2 nm, do = 9 nm, As-t =125 /reV and zero mag¬ 
netic field. The plotted energy levels at the working point 
(F = 2.73 mV/nm), denoted by a dashed line, correspond, 
from bottom to top, to the modified singlet state | S)\ the 
degenerate triplet states |T_), |To) and \T+) and two linear 
combinations of the two doubly occupied states | Rf, 0} and 
|0,4t) with small (5.4% and 0.05% at the working point) sin¬ 
glet state admixtures. The lowest/highest energy levels were 
shifted down/up by 0.5 meV for the energy gaps to become 
visible in the plot. These energy gaps are, from top to bot¬ 
tom, the anticrossing splitting of the doubly occupied states 
which equals 35 /reV and the singlet-triplet splitting which is 
35/reV at the center of the splitting (F = 3.71 mV/nm) and 
125/xeV at the working point. 


the creation of an exciton at the bottom dot by the laser 
field. The 4 trion states have the same single-particle and 
Coulomb interaction terms in their energies. Moreover, 
we calculated for a wide range of relevant experimental 
values that the electron-hole spin exchange interaction 
is negligible, since the unpaired electron and hole reside 
in different dots. We also found that hole tunneling is 
too weak to have a notable effect on the energies due to 
the highly localized hole wave functions. We therefore 
conclude the trion states are degenerate. 

With the energies of the spin states and the trion states 
determined, we obtain the selection rules and plot the 
energy level diagram for the system in the Faraday ge¬ 
ometry and zero magnetic field in Fig. ^|. The allowed 
transitions are with circularly polarized light with po¬ 
larizations a±. We notice the \T±) states are each iso¬ 
lated from the rest of the spin states in the sense that no 
Raman transition can be induced between each of these 
states and another spin state. This result, which persists 
when a magnetic field in the 2 direction is applied, pre¬ 
cludes complete optical control of the two-electron spin 
state in this geometry and, consequently ,we decide to 
consider the Voigt geometry with a magnetic field ap¬ 
plied in the +x direction. 

The application of a magnetic field in the +x direction 


m,t) 


IttU) 


m,t) 


m,i) 



FIG. 3. (color online) Energy level diagram of the quan¬ 
tum dot molecule in the Faraday geometry with zero magnetic 
field. The states include the modified singlet state |S); the de¬ 
generate triplet states |To), |T+) and |T_); and the optically 
excited degenerate trion states | 41iht)> I 4-tlh 4) i I 4t4it) 
and | htfli-l) with energy Et r ■ The kinetic exchange interac¬ 
tion shifts the modified singlet energy level below that of the 
triplet states with the splitting magnitude being A s-t- The 
selection rules are shown in the diagram with the red/blue 
arrows denoting u± polarization. 


Zeeman-splits the triplet states in energy to the states 
\T-) X , \T 0 ) X and \T+) X as shown in Fig. 0 while the 
modified singlet state, having no spin projection in the x 
direction, is not shifted in energy. The new triplet states 
are given by 


where 


\ T -)x = 

(18) 

|T 0 ), = 2- 1 / 2 (|+, — ) + | — , +> ), 

(19) 

\T+) X = |+, +}, 

(20) 

l±>= 2- 1/2 (| t)±U» 

(21) 


are the eigenstates of the spin operator in x direction. 
The trion states are also Zeeman-split in energy, with 
the new states being |f + ,±) and |t_,±) and where the 
states |t±) are defined by 

|f±> = 2- 1/2 (| 4tlT) ± I m)- (22) 

We plot the trion states Zeeman splittings in Fig. [|. With 
the energy levels determined, we write the selection rules 
for the Voigt geometry. We find the allowed transitions 
have linear polarizations in either the x or y directions, 
and that the set of transitions, shown in Fig. [|, makes Ra¬ 
man transitions possible between all pairs of spin states 
apart from the pair of |T_) x and \T+) X , a point that will 
be important later when we tackle the problem of com¬ 
plete quantum control of the two-electron spin state. 
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FIG. 4. (color online) Energy level diagram of the quantum 
dot molecule in the Voigt geometry (growth direction is z, 
and magnetic field is in +* direction). The energy of the 
modified singlet state |>S) is unaffected by the magnetic field 
and remains split from |7o) x by As-t, while the triplet states 
are Zeeman-split in energy to \T + ) X , |To) x and | T-) x , with the 
splitting being A = |gr e |/rsI3, g e the electron p-factor, g ,b the 
Bohr magnetron and B the magnitude of the magnetic field. 
The degeneracy of the optically excited states is also removed 
by the magnetic field and the states are Zeeman-split as shown 
in the figure, where A^ = |< 7 /i|ps-B and gh is the hole g-factor. 
The selection rules are plotted in the diagram with red arrows 
corresponding to V ( 7 t v ) polarization, and blue arrow arrows 
to H (tv x ) polarization. 


III. UNIVERSAL COMPUTATION WITH 
TWO-LEVEL OPERATIONS 

The spin states in the quantum dot molecule system 
in the Voigt geometry form a basis for a two-qubit com¬ 
putational state space. This basis is |±,±), where the 
|±) states were defined in Eq. (|2l|), and where the ad¬ 
mixtures of the doubly-occupied states were omitted for 
brevity. The scheme for universal computation we pro¬ 
pose provides a realization of a given arbitrary unitary 
operation in the computational basis through the appli¬ 


cation of a set of two-level unitary operations between 
the spin states. Each of the two-level unitary operations 
operates on two spin levels and is realized through two 
laser pulses with the same Raman detuning S that pump 
the transitions from the two spin states to a common 
trion state. 

Suppose we want to realize a given a unitary operation 
U in the computational basis as a product of two-level 
unitary operations. A general proper unitary two-level 
operation in the space of the states |i) and |j) ( i,j = 
1,...,4) in the spin state basis , |S), \T 0 ) X , \T + ) X and 
\T-) X , is given by 

Q 

Rij(6,n) = exp(-i-n ■ <x), (23) 

which rotates the pseudo-spin vector in the subspace of 
the states |*) and | j) by an angle 9 about the axis n, 
with <j = ( <J x ,cty , a z ) acting in the subspace of these two 
levels. Going back to the unitary operation U, we first 
write its matrix relative to the spin state basis. Then, 
Ref. |3^ shows how to decompose an arbitrary n-qubit 
unitary operation to a product of up to ( 2 9 ) two-level 
unitary operations. However, we prefer the operations in 
the product to be proper unitary, as such-aperations are 
more easily implemented experimentally, E-a and we give 
in Appendix |Aj the procedure for decomposing U to a 
product of up to 6 proper unitary two-level operators in 
the form of Eq. ( P3| ) and an overall phase factor, which 
can be omitted. 

From the energy level diagram in Fig. [| it can been seen 
that apart from R34 each of the 6 two-level unitary oper¬ 
ators that may appear in the decomposition of U, namely 
R12, R13, Ru, R23 , R24 and .R 34 , may be implemented 
with a single Raman transition. A graphical illustration 
of the situation appears in the graph in Fig. ^], wherein 
each vertex represents a level and each edge a possible 
two-level unitary operation. The operation R34 may be 
implemented through 3 Raman transitions as shown by 
the following argument. If we represent a 7 r rotation op¬ 
eration between the states |«) and |j) as 

Pij = Rij(n,y), (24) 

then we have 

Rik = PjkRij Pjk ’ ( 25 ) 

Rkj = PikRijP}k' (26) 

With the problem of R34 solved, we consider the fact 
that it may be advantageous to use of only some of the 
two-level unitary operations, as the realizations of those 
using Raman transitions may be faster or have higher fi¬ 
delity relative to the other operations. We therefore ask 
what universal subsets of the 5 operations can be chosen, 
i.e. subsets with which any arbitrary unitary operation 
can be implemented. The answer to this question lies in 
the graph interpretation of Eqs. (p5[j2^). These equations 
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FIG. 5. Graph representation of the spin states and the 
possible two-level operations between them. The graph ver¬ 
tices represent levels and its edges represent possible two-level 
unitary operations. The vertices 1, 2, 3 and 4 are associated 
with the states 15), |To}a;, | T + } x and |T_) X , respectively. 


show that an edge may be realized by 3 Raman transi¬ 
tions if a two-edge path connects its vertices. Extending 


this result by induction, we find that an edge can be re¬ 
alized by 2n — 1 Raman transitions if a path of n edges 
connects its vertices. An important result follows. For 
a subset of the graph in Fig. |] to be universal, it must 
include all vertices and be connected. When we define n- 
edges as edges representing ± 7 r rotations, the graph for 
the minimal universal set of operations can be found and 
its form is given in App. |b|. 

IV. PULSE SHAPING 

The two-level unitary operations in our system are re¬ 
alized through Raman transitions. For such a transi¬ 
tion to implement the quantum dot molecule is il¬ 
luminated by two phase-locked laser pulses propagating 
in the 2 direction. The two pulses are linearly polar¬ 
ized, following the selection rules in Fig. |[ in either the 
vertical ( V ) or horizontal ( H) directions, have a com¬ 
mon Raman detuning S and pump the two transitions 
from the spin states |?) and | j) to a common trion state. 
The model Hamiltonian in the basis of |S), \Tq) x , \T + ) x , 
\T-) X , |f+,+), 1 1 + ,~), |f_,+), and |t_,-), in the rotat¬ 
ing wave approximation (RWA) is 


/- 

-As_t 

0 

0 

0 

X*H 

—Xv 

* 

Xv 

% 

-Xh 

\ 


0 

0 

0 

0 

Xv 

~X*H 

X*H 

—Xv 



0 

0 

-A 

0 

V2 X *v 

0 

V2 X *h 

0 



0 

0 

0 

A 

0 

V2 X *h 

0 

V2 X * v 



XH 

Xf 

V%xv 

0 

Etr + {Ah — A)/2 

0 

0 

0 



-Xv 

~XH 

0 

V%XH 

0 

Etr + {Ah + A)/2 

0 

0 



Xv 

XH 

V^XH 

0 

0 

0 Etr 

- {Ah + A)/2 

0 


V 

- XH 

-Xv 

0 

V?XV 

0 

0 

0 Etr 

— {Ah — 

A)/2 / 


(27) 


where E tr is the trion state energy shown in Fig. [i], we 
have set U = 1, A = \g e \nBB is the electron Zeeman 
splitting, A h = \gh\fJ-BB is the hole Zeeman splitting, 
and where /i n is the Bohr magnetron, B is the magnetic 
field magnitude, and g e and gh are the electron and hole 


g- factors, respectively. The expressions xv(t) = fV(t)/2 
and X-ff(i) = fiff(t)/2 in Eq. (^?|) are, respectively, halves 
of the time-dependent Rabi frequency for the transitions 
\S) -A |f_,+) and |Xo)x ~ 

In the interaction representation the Hamiltonian in 
Eq. (^?]) is recast as 
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(28) 


where A,; ? = Ei — Ej , and where Ei is the energy of 
the state |«). Eqs. (|27]) and (p8|) were derived with | S) 


approximated as the singlet state | S). When the full 
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expression for IS} is used, the Hamiltonian matrices are 
the same, apart from an extra numerical factor in the 
off-diagonal matrix elements involving | S). This factor is 
the coefficient of |S) in |S), which for the experimental 
values in Ref. |2^ is 0.973. 

As seen from the energy level diagram in Fig. [|, a 
Raman transition from one spin state to another may 
suffer from unintended dynamics due to the presence of 
close resonances. Rather than increase the pulses lengths, 
which will make the operation slower and magnify the 
effects of dephasing, we copejsjith the unintended dy¬ 
namics through pulse shaping, eSD-j whereby we vary the 
parameters of the pulses to optimize the operation. 

We describe the explicit form of \v(t) and XH(t) with 
pulse shaping for the V-H, H-H and V-V cases, with the 


case names denoting the polarizations of the Raman laser 
pulses involved. In the V-H case we have 

Xv it) = xo(t) + X 2 (t), (29) 

XH(t) =Xi(t- Af) + x 3 (f- At), (30) 

where 

Xj(t) = Xje- {t/Tj)2 e~^ t+i ^, (31) 

Xj{t) = Xj(t)e iv 'i lnc ° sh (‘/T 3 ) ; (32) 

and where \j(t) (j = 0 , 1 ) is a main pulse with ampli¬ 
tude Xji pulse width r ? . central frequency uij and con¬ 
stant phase ( pj. The shift in time between the two main 
pulses is given by At. The expression Xj{t) ( j = 2,3) 


in Eq. (g2j) is a helping pulse designed to alleviate the 
effects of unintended dynamics. This pulse has an ad¬ 
ditional chirping termJiU which sweeps the instantaneous 
frequency in time through a range of frequencies centered 
about u>j. 

The pulse shapes for the other two cases involve only 
one polarization type. In the H-H case we have 


Xv(t) = 0, (33) 

XH(t) = Xo(t)+X2(t) + 

Xi(i- Af) + x 3 (i- At), (34) 

and in the V-V case we have 

Xv (t) = Xo(t)+X 2 {t) + 

Xi{t-At) + X3(t-At), (35) 

XH(t) = 0. (36) 


The pulse parameters should be varied so as to max¬ 
imize the fidelity of the operation relative to the ideal 
Rij. As this fidelity depends on the input state, the ex¬ 
pected value of the fidelity over all possible input states 
is required. Let the states in the basis | S), \Tq) x , \T + ) x , 
\T-) X , |£+,+), |i+,-), |i-,+), and |f_,-) be denoted 
by | j) ( j = 1,...,8). For unitary evolution we find in 
App. 0 the expected fidelity is 


where 


hj = (i\&U id \j), (38) 

and where U and Uid are, respectively, the actual and 
ideal evolution operators. The expected fidelity in 
Eq. (p7|) is the yardstick by which we decide whether 
one set of pulse parameters is better than another set of 
parameters for the implementation of a given two-level 
unitary operation Rij. 

We now incorporate the effects of light hole mixing in 
our analysis. The hole states | -O') and | 1}) in Eq. (|2^), 
denoted as \Hf), and heretofore taken as bare heavy hole 
states, are better approximated through the Luttinger 
HamiltonianliZlua as a superposition of heavy hole states 
and light hole states of the form 

I H t) = cos0 m ||,±|) - sin0 m e =F ’*™|!,=F^). (39) 

In Eq. (H|) 6 m and </> m are the mixing angles, ||, ±|) are 
the heavy hole states, ||,±^) the light hole states, and 
the quantization axis is taken as the growth direction (z 
direction). 

When the new hole states in Eq. ([hi]) are substituted in 
the trion states and the Hamiltonian matrix is calculated 
following an analogous derivation in the single dot case 
in Ref. pi we find results identical with the ones obtained 
in this reference, namely that adjusting the polarizations 
of the lasers from V and H to V' = 2- 1 / 2 (cr + + e v +cr_) 
and H' = 2 ~ 1/2 (cr + - e~^-(T_), with 


e i M± = VScosOm ±smd m e ±i(l,m 
v/3 cos 9m ± sin 9 m e ’ 

the Hamiltonian has the same form as in the case with 
no light hole mixing [Eq. (^ 8 |)] apart from replacing xv 
and Xh by 

1 + 2 cos 2 6 m 

3 cos 6 m + sin 9 m ’ 

1 + 2 cos 2 9 m 

Xh Xh ^ cQg ^ g j n ’ 

where Xv and xh' are the analogous quantities to Xv 
and xh for light with V' and H' polarizations, respec¬ 
tively, and the effects of V' and H' being non-orthogonal 
are neglected. 

The identity of the last results with the single dot case 
in Ref. H is not surprising, as only the bottom quantum 
dot is optically excited. To summarize, the effects of light 
holes mixing can be circumvented by adjusting the polar¬ 
ization of the laser pulses and their Rabi frequencies as 
noted above. With these adjustments, the Hamiltonian 
in Eq. (p8|) can still be used. 


(41) 

(42) 


V. DECAY AND DECOHERENCE 


^ = + I 2 ), 


*=1 3=1 


(37) 


When decay and decoherence are taken into account 
the state of the system may be described by a density 





matrix p and its non-unitary evolution by a quantum 
master equation in the Lindblad form 

p=-i[H,p] + 'Y^£i[p]> ( 43 ) 

i 

where H is given in Eq. (|2^) , and where the sum is over all 
trion state relaxation channels, each of which is described 
by a Lindblad superoperator 

C t [p\ = D iP D\ - \D\D iP - \pD\D (44) 


In taking the relaxation channels as separate and inco¬ 
herent, we neglected spontaneously generated coherence 
(SGC), since we assumed the Zeeman splittings are large 
enough to satisfy |Ayfc| F (j,k = 5,...,8; j ^ k), 
wherp-|E-. is the total relaxation rate of a given trion 
level,and that the level splittings A, A h and A s~t 
are_larger or of the same order as the Rabi frequency 
Q.El The effects of pure dephasing are not included in 
Eq. (|43|) , as this process, which is generated by the nu¬ 
clear spins, has a rate much kmipr than the total trion 
dephasing rate in this system.0L3 

The 12 jump operators, D i: in Eq. (43), each corre¬ 
sponding to a possible relaxation channel, are 


Di = Vm\S)(t+,+\, 

(45) 

D2 = y/T/i\T 0 ) x (t+,+\, 

(46) 

D 3 = y/f/2\T+) x (t+,+\, 

(47) 

D 4 = y/TJi\S)(t+,-\, 

(48) 

D 5 = ^/Tj4\T 0 ) x (t + ,-\, 

(49) 

D 6 = y/T/2\T-) x {t+,-\, 

(50) 

D r = y/fJi\S)(t-,+\, 

(51) 

D 8 = y/T/i\T 0 ) x {t-,+\, 

(52) 

D 9 = y/fj2\T+) x (t^+l 

(53) 

D w = ^/Tj4\S)(t-,-\, 

(54) 

D n = y/r/i\To) x {t-,-\, 

(55) 

D 12 = y/T/2\T-) x {t-,-\. 

(56) 


The terms in the square roots in Eqs. ( ph||56| ) are the 
relaxation rates for the corresponding decay channels. 
These rates were determined by splitting the total decay 
rate, the trion state linewidth T, in the ratio of absolute 
square values of the dipole matrix elements for the al¬ 
lowed transitions from the this state to the spin states, 
as these values are proportional to the spontaneous decay 
rates of the corresponding channels. 

With spontaneous decay and decoherence included in 
the model, the expression for the expected fidelity of an 
operation relative an ideal operation in Eq. ( |37| ) should 
be generalized. The derivation in App. [d] gives this gen¬ 
eralized expected fidelity as 


T = 25 E E 01 U&pWUuli) + ml d p^U id \j) 


*=1 3 =1 


( 57 ) 


where pbi) is the density matrix resulting from the non- 
unitary evolution of an initial density matrix |i)(jj. In 
the limit of unitary evolution (T = 0), pW) = U\i)(j\U^ 
and Eq. (|57j) reduces to Eq. (|37|) . 


VI. THE TWO-QUBIT QUANTUM FOURIER 
TRANSFORM 


As a case in point of a non-trivial two-qubit transfor¬ 
mation realized using the methods described above, we 
choose the two-qubit quantum Fourier transformation. 
In the computational basis |±, ±) this transformation is 


Uqft = 


1 



(58) 


Rewriting Eq. (|5^) relative to the physical basis | S), 

| T 0 ) x , | T + ) Xl and \ TJ) X as Uqft , and applying the algo¬ 
rithm in App. |a] to decompose the result to a product of 
two-level proper unitary operations, we find 

Uqft =e lQ i? 32 ( 032 )^32)^34(^34?^34)^12(^12)^12) x 

7 ^ 14 (^ 14 , ^ 14 ) 7 ^ 24 (^ 24 ? ^ 24)5 (59) 


where Rij(9,n) is defined in Eq. (p23|) , and where a, 6 ij 
and riij were numerically determined. 

We detail the optimization process of the pulse train. 
First, we choose the physical parameters for the system 
as the ones in Ref. [2^, namely hs = 2.6 nm, h F = 3.2 
nm, do = 9 nm (see Fig. [l]) and Ag-T = 125 peV. 
The electron and hole p-factors are taken from Ref. |] as 
g e = —0.48 and gh = 0.31, and the trion states linewidth 
r is obtained from Ref. as T = 1.2 peV. Next, for each 
of the operations R 12 , R 13 , R 14 , R 23 and R 24 we choose 
the trion state that will be used to realize the Raman 
transition. We choose this state so as to reduce the unin¬ 
tended dynamics, i.e. such that the resonance frequencies 
of each of the transitions from the two spin states to this 
trion state is as far as possible from other allowed tran¬ 
sitions with the same polarization. We then choose a 
value for the magnetic held B that maximizes the min¬ 
imum absolute frequency difference among all pairs of 
allowed transitions with the same polarization. For mag¬ 
netic held values lower than 10 T, this value is found to 
be B = 8.99 T, and the minimum absolute frequency 
difference is 36.5 peV. 

We then turn to the optimization of the individual 
pulses. Since each Ri :] in Eq. (59) (apart from R 34 ) can be 
implemented directly as given or indirectly via Eqs. ( p5| - 
pMi| ) , we consider all possible indirect implementations and 
optimize the 35 operations Rij(O 32 , 7432 ), Rijifi 34 , 7434 ), 
Rij{0i2,ni 2 ), Rij (014,flu), Rij^nu), Rij and Pk , 
where (i,j) = (1,2), (1,3), (1,4), (2,3), (2,4). In each op¬ 
timization we realize the operation with two main pulses, 
each of which having the form of Eq. (Bh), and maximize 
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the expected fidelity in Eq. ( p7| ) through grid search in 
pulse parameter space followed by gradient descent. In 
the optimization process we keep in mind the need for 
short pulses that will ensure fast operation. Examining 
the results of the optimizations, we conclude that direct 
implementation is always preferable when possible, with 
the only case where indirect implementation is chosen 
being 

^ 34 (^ 34 ,^ 34 ) = -Pl3-Rl4(#34, n34)Pi 3 , (60) 

which cannot be directly implemented due to the selec¬ 
tion rules. 

With the pulse train determined to consist of 7 oper¬ 
ations, we add two helping pulses, each in the form of 
Eq. (p^), to each operation. The parameters of the first 
helping pulse are varied so as to maximize the fidelity 
in Eq. (£7|) and the process is repeated for the second 
helping pulse. Following the optimization of the helping 
pulses, we turn to the optimization of the time spacings 
between the 7 operations. The fidelity we maximize is 
now the fidelity of the entire QFT with decay and de¬ 
coherence taken into account as given in Eq. (]5?|). The 
resulting optimized pulse trains are calculated to have an 
expected fidelity of 84.9% and a duration of 453 ps when 
helping pulses are not employed, and an expected fidelity 
of 88.1% and a duration of 414 ps when helping pulses 
are used. The main source of the infidelity is unintended 
dynamics, that is excitation of levels other than the ones 
intended. This result is revealed comparing the fideli¬ 
ties in Eq. (|37|) and Eq. (|t|) for each of the individual 
operations. We find the decrease in fidelity from unity 
due to decoherence and decay is, on average, 0.8% per 
operation, while the decrease due to unintended dynam¬ 
ics is 2.6% per operation. This result could have been 
expected as the minimum absolute frequency difference 
of 36.5 /reV is always comparable to or smaller than the 
spectral width of the pulses. 


VII. CONCLUSIONS 

In this work, we demonstrated a scheme for complete 
quantum control of a system of two singly-charged QDs 
coupled by coherent electron tunneling. We derived the 
energy level diagram for the system and chose the two- 
level operations realized by the 5 possible Raman transi¬ 
tions through the trion states as the universal set of gates. 
In contrast with the commonly used universal computa¬ 
tion scheme of CNOT gates and single-qubit gatesj23 we 
used only two-qubit gates, each realizable through a few 
laser pulses. Our choice is advantageous as single-qubit 
gates, which can be described as two simultaneous two- 
qubit gates, r-tcnd to suffer from unintended dynamics in 
this system.L3 We demonstrated our scheme by describ¬ 
ing the realization of the two-qubit QFT in the QDM 
system with a reasonably high fidelity of 88.1%. We in¬ 
corporated light hole mixing, decay and decoherence in 
our analysis. 


Our results indicate that pulse shaping can substan¬ 
tially increase the fidelity of quantum operations in opti¬ 
cally active QDs, and that the QDM is a promising plat¬ 
form for multi-qubit quantum computation. The main 
limitations for increasing the overall fidelity of a given 
two-qubit operation with this scheme are unintended dy¬ 
namics, decay and decoherence. The first can be dealt 
with through improved pulse shaping. Due to limited 
computing power, we chose our pulse search algorithm 
to optimize each pulse separately. A better optimization 
may be achieved through simultaneous optimization of all 
4 pulses of a given operation. Ameliorating the effects of 
decoherence and decay can also be achieved by improved 
pulse search methods. In our search we maximized the 
fidelity expression in Eq. (|37]) that does not take decay 
and decoherence into account, rather than the expression 
that takes these effects into account in Eq. (^7|). Using 
the latter expression, which takes much longer to evalu¬ 
ate, will increase the fidelity of the overall operation. 

The experimental realization of universal computation 
in QDMs will be a major step towards building a QD- 
based optically-controlled quantum computer. Future re¬ 
search may look into methods for robust and effective 
entanglement generation between distant QDMs, possi¬ 
bly wilh_.multi-photon light, as was suggested for single 
QDs.Ela Quantum computations with more than two 
qubits may be realized in a quantum network of entan¬ 
gled QDMs through unitary operations on each of the. 
QDMs, possibly employing the theory of cluster states.E3 


Appendix A: Decomposition of a Unitary Operation 
to a Product of Two-Level Proper Unitary 
Operations 


In this appendix we show a procedure to decompose 
a given unitary operation U to a product of an overall 
phase factor e la and two-level proper unitary operations. 
The procedure is given for the case of a 3-state basis, but 
can be easily extended to an V-state basis with N > 3. 

The first step in the algorithm is to factor out \U\ = e la 
and proceed with U' = e~ la ^ 3 U, which is proper unitary. 
Let U' be given by 


U' = 




(Al) 


If b = 0, we take U\ as the identity matrix and proceed 
to the next step. Otherwise we take U\ as the proper 
unitary matrix 


Ui = 


1 

vw+w 


a* b* 0\ 
—b a 0 j 
0 0 1 / 


(A2) 


and have 


U X U' 


a' d! g' \ 
0 e' h'\ 

c' S’ f j 


(A3) 
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Then, if d = 0, we take C /2 as the identity matrix and 
proceed to the next step, while if d 7 / 0 , we define 


U 2 


1 


a'* 0 d* \ 

0 10 , 

-d 0 a' ) 


(A4) 


which is proper unitary, and have 

/! 0 0 \ 

U 2 UiU' = 0 e" h" = U 3 . (A5) 

V 0 f" j" ) 

Hence, the sought decomposition is 

U = \U\ 1/3 Ud 1 U2 1 U 3 . (A 6 ) 


Appendix B: The Minimal Universal Set of 
Operations in the Quantum Dot Molecule 


In this appendix, we prove the minimal universal set of 
two-level unitary operations in the quantum dot molecule 
system is represented by the graph in Fig. In the last 
paragraph of Sec. Ill it was shown that a graph repre¬ 
senting a universal set of operations must be connected 
and contain all vertices. Hence, the minimum number of 
edges is 3. The edges cannot all be it edges since then the 
number of possible two-qubit unitary operations attain¬ 
able with a finite number of such two-level operations is 
finite. We conclude the set of operations represented by 
the graph in Fig. is minimal. 

This set of operations is universal, since each of the 
6 two-level unitary operations that can appear in the 
decomposition of an arbitrary unitary operation U can 
be written i n term s of the 3 operations in the set through 
, namely 


use of Eqs. (|25 


with Rij a 
spectively. 


Rl2 

= -R 12 , 


(Bl) 

Rl3 

= P 23 -R 12 -P 23 ! 


(B2) 

Rl4 

= P21P14R21P14 

P t 
-*21 ? 

(B3) 

R23 

= P12P23R12P23 

p t 

-*12 5 

(B4) 

R24 

= P14R21P14, 


(B5) 

R34 

= P14P23R21P23 

p t 

^14’ 

(B 6 ) 

Pij defined in Eq. (|23|) and Eq. (|24| 

), re- 


Appendix C: The Expected Fidelity for Unitary 
Evolution 

In this appendix, we derive the expression for the ex¬ 
pected fidelity of an operation relative to an ideal oper¬ 
ation when state evolution is unitary. The input state 
l^in) is given by 

4 

li’in) =^2bj\j), (Cl) 

i=i 



FIG. 6 . Graph representation of the minimal subset of the 
possible two-level operations that can still achieve universal 
computation. The graph vertices represent levels and its edges 
represent possible two-level unitary operations. A 7 r-edge rep¬ 
resents a ± 7 r rotation operation. The vertices 1, 2, 3 and 4 
are associated with the states \S), \To) x , \T+)x and | T-) x , 
respectively. 


while the output state is 

8 

\i>out) = U\tpin) = ^2 C J'U)> ( C2 ) 

.7=1 

with U being the evolution operator. The ideal output 
state is given by 


4 

I Ad,) = U id \ip in ) = ^2 d i L?)> ( C3 ) 

. 7=1 

where Uid is the ideal evolution operator. The fidelity 
of the output state relative to the ideal output state is 
calculated as 


P = {'4>id\Tr T {\ip out XV* out 1)1 Ipid), 


(C4) 


where we trace over the trion states. Plugging Eqs. ( |C2| ) 
and (CS) in Eq. (C4), we find 


T = 


12 c *d 




3 =1 


(C5) 


Averaging Eq. ( |C5| ) over all possible input states as in 
Ref. 44, the expected fidelity of the operation is obtained 
in the form of Eq. f[i7|). 


Appendix D: The Expected Fidelity for 
Non-Unitary Evolution 

In this appendix, we derive the expected fidelity of 
an operation relative to an ideal operation when state 
evolution is non-unitary. Let the initial state be a pure 
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state |i pi n ). This state evolves to the density matrix p 
with the evolution written in terms of effect operators 
Ak as 


JVi 


P = '^2 Ak\lf>in){lpin\Al. 


(Dl) 


fc =1 


The fidelity of p in Eq. (Dl) relative to an ideally evolved 
state | il> id ) = U id \ip in ) is given by 


JVi 

T = J2\(^n\A{Ui d \^ in )\ 2 . (D2) 

k =1 

When each term in Eq. (|D^) is averaged over all possible 
input states, we find, using Eqs. (flg-f^l). the expected 
fidelity is 


N i 4 


k—1 2—1 j =1 




0)|2 


(D3) 


where 1^ = {i\A^ k Ui d \j). Changing the order of summa¬ 
tion in Eq. ( |P3| ) and using Eq. (Dl) for the time evolu¬ 
tion, we obtain Eq. (|57|). 
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